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We have calculated a discrete spectrum and found an exact analytical solution in the form of 
Meixner polynomials for the wave function of a thin gravitating shell in the Reissner-Nordstrom 
geometry. We show that there is no extreme state in the quantum spectrum of the gravitating shell, 
as in the case of extreme black hole. 
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In the absence of a quantum theory of gravity, the 
semiclassical approximation (gravitons, the Hawking ef- 
fect, the primordial perturbation spectrum in cosmology) 
and various semi-qualitative ''plausible" models (for the 
black hole mass spectrum) [l[ Q are used as a convenient 
method for describing quantum effects against the classi- 
cal background of general relativity. One of the models, 
in which it is relatively easy to obtain the black hole mass 
spectrum, is the model of a thin gravitating shell Q . 

There exists a simple method for obtaining the mass 
spectrum in the formalism of thin shells. It is based on 
the natural assumption that the mass of a gravitating 
system (e. g., a black hole) at spatial infinity, m out , is 
the Hamiltonian of the system, because m ou t is the to- 
tal energy of the entire system that is conserved during 
the dynamical evolution of the shell. Given the Hamilto- 
nian, the wave equation can then be easily written and 
the corresponding mass spectrum can be found from its 
solution. 

All properties of a Schwarzschild black hole in general 



relativity are known to be completely determined by one 
parameter, namely, the mass at spatial infinity. However, 
if we consider the Carter-Penrose diagram for an eternal 
Schwarzschild black hole (see the figure 1), which de- 
scribes the global geometry of the manifold in question, 
then we will see that there exist two spatial infinities, 
the so-called regions i?+ and Therefore, it can be 
assumed that the mass spectrum in the Schwarzschild 
metric of an eternal black hole can depend on two quan- 
tum numbers 

Consider the dynamical equation that describes the 
evolution of a thin gravitating shell in the Reissner- 
Nordstrom metric. Thi s eq uation has a well known form 
(see, e. g., 0,Q and O-Ql): 
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= 47rpju(p), (1) 

P P 

m ou t, Qin and Q out are, re- 
spectively, the black hole mass and charge inside and out 
side the shell; p = p(r) is the shell radius relative to an 
observer on the shell; and the dot over the function de- 
notes a derivative with respect to the proper time r of 
this observer. We will consider only a dust shell, then 
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the function p(p) = A/p 2 , where A > is the constant 
of integration [6|. For the subsequent analysis, it is con- 
venient to define a quantity M — 4nA that, as can be 
shown, is the total mass of the shell. As is easy to show 
from Eq. ([!]), the conditions for the signs of <7i njOU t are 
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Everywhere below, we will consider the case where 
m out > m in , M 2 > Ql ut - Qf n , and cr in = 1 in the entire 
spacetime region. At the same time, the sign of cr ou t can 
have any value. Squaring Eq. ([1} yields an expression of 
the form 
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The total mass of the entire system m out , which is con- 
served with time, appears on the left-hand sid e of this 
expression. Following the earlier papers 
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13j, we will 



define this mass as the Hamiltonian of our system H (see 
also where a slightly different model was consid- 
ered). Let us make the change of variable x = Mp in 
Eq. (4} (see 
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The Hamiltonian will then take 



the form 



H = a in ^x 2 +M 2 [l 

+m in + M^ t ~' : '" 
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Given Hamiltonian ([5]), we calculate the Lagrangian of 
the system under consideration 



• In 



+ Jx 2 + M 2 1 - 



2m iB M , Q 2 n M 2 



*Jx 2 + M 2 (l- 2 -^ + 9lMl 

° ut + ^ in <7 in x\nM-m in , (6) 



+M 



M 2 Q 2 out + Ql 



2x 




Figure 1: Carter-Penrose diagram for an eternal 
Schwarzschild black hole. 



and then the canonical momentum 



p = a in \n 
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Let us now write Hamiltonian (JS|) in terms of the canon- 
ical momentum ([7]). For this purpose, we will find the 
variable x from Eq. ([7]) and substitute it into Hamilto- 
nian ([S]). As a result, the dependence of the Hamiltonian 
on the coordinate x and canonical momentum p will take 
the form 



H = min - M 



M 2 - Q 2 out 



2x 



(8) 
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Let us write the wave equation H<j)(x) — m ou t0(a;), using 
a standard commutation equality of the form [p, x] = —i 
and the identity 



exp (^o 4>{x) = 4>{x + x ). 

As a result, the wave equation will be written as 
2m ln M , Qf a M 2 



(9) 



(x-i) + 1- 



x + i) 



m 2 - n 2 4- o 2 

* aat + Vm <l>(x)=2E<t>(x), (10) 



where E = (m out — m in )/M. The wave equation (|10|) 
differs from the standard Schrodinger equation in that 
this is not a differential equation but a difference one. 
This is because the quantization is performed over the 
proper time of an observer on the shell and not over the 
time of an observer at infinity. If we represent the expo- 
nential as a series, then we will obtain an infinite order 
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differential equation. Consequently, this equation should 
be supplemented by infinite boundary conditions. These 
were found in [3] and, in our case, are 



0^(0) = 0, 1 = 0,1.. 



(11) 



It should also be required that the wave functions do not 
diverge at spatial infinity. 

The solution of the wave equation can be expressed in 
terms of Meixner polynomials (see, e. g., [ljj]) that satisfy 
the equation 

a(x)[f(x + l)-2f(x) + f(x-l)} 
+r(x)[f(x + l)-f(x)}+Xf(x)=0, (12) 

where a(x) — x, a(x) + t(x) — ^(7 + x), and /1 and 7 
are some constants. We will solve the wave equation (fpQ|) 
in an explicit form by reducing it to an equation of form 
(fT2j) . To begin with, let us make the coordinate trans- 
formation x — > —ix which means the rotation through 
— 7r/2 in the complex plane. In this case, Eq. (flO| will 
take the form 

2m in Mj Q?M 2 " 



<f>(x + l) + 1 ^^U(x-l) 

\ X x A ) 

_ ^-Q^ + Ql ^ = 2E4>{x). (13) 

X 

First, consider the simple case where m m = Q m = 0. 
The solution of Eq. (IT31 can then be expressed in terms 
of Meixner polynomials and is 

(3 2x T{x) 
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where 



13 = E + ^E 2 - 1, Af(x) = f{x + 1) - f(x), 
C{x) = C{x + 1), (15) 

T(x) is the gamma function, and C{x) is a periodic func- 
tion with a period of 1. Let us expand the periodic func- 
tion C(x) into a Fourier series: 



C( x ) — Ck e~xjp(2irikx) . 



(16) 



Exactly the same factor appears in the solution of the 
wave equation if we pass from the coordinate representa- 
tion to the momentum one H] • The coefficients 
can be found from the boundary conditions to the wave 
equation (fTT|) . It is easy to show that the wave functions 
4> n {%) are orthogonal in the following sense. If we set 
Xi = x and Xi + \ = Xi + 1, then we can write the follow- 
ing sum using the properties of Meixner polynomials: 
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(17) 



where the weight function p[x) = l/[xC 2 (s)] and 
2 T(n)T(n + 1) 

It follows from Eq. (jTTJ) that the wave functions 4>n(x) 
at n 7^ m are orthogonal, provided that < /3 < 1 [1 51 ] . 
These wave functions are easy to modify in order that 
they also be orthonormal. 

The discrete mass spectrum E n corresponding to the 
wave functions 4>n(x) in the form of Meixner polynomials 
15j satisfies the simple equation 



i(M 2 - Q 2 out ) = 2ny/E 2 - 1. 



(18) 



As a result, for the sought for mass spectrum we find 
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The hydrogen-like mass spectrum found generalizes the 
result obtained previously in [ill ] and coincides with the 



result of 



13, 



1G, 



17|, where it was obtained by a dif- 
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ferent method. A discrete mass spectrum takes place if 
E n < 1 . If, alternatively, E n > 1 , then a continuous mass 
spectrum will take place [ijj]. Note that, in accordance 
with Eq. (|T9"| for the mass spectrum, the total gravitating 
mass of the metric m ou t turns out to be lower than the 
classical "bare" value of M due to the quantum correc- 
tions, as it must be for a gravitationally bound system. 

When (M 2 — Qout) 2 /4 > 1, we can write a condition 
for the applicability of the semiclassical approximation. 
Indeed, in this case, there exists a minimum value of 
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the quantum number, n m i n = \(M 2 — Q 2 ut )/2], where 
" [] " denotes the integer part, for which the semi-classical 
approximation will definitely hold at n » n m i n . In the 
case, where the opposite inequality (M 2 — Q 2 ut ) 2 /4 < 1 
holds, the semiclassical approximation is applicable at 
any n. 

To write the solution of the original equation (fTU)) . the 
inverse change of variable x — > ix should be made. As an 
example, let us write out the first several polynomials: 



(x) n=1 = ix(3 l 
(x) n=2 = ix/3 l 
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Ck exp(— 2irkx), 



1 



> Cfeexp(— 2irkx). 
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When the requirement that the wave functions be finite 
at infinity is fulfilled, the constants = will be zero 
for all k < 0. The remaining constants of integration 
can be easily found from the boundary conditions at the 
coordinate origin. This was done in [lH for n = 
1, and we will not give them here. Making the change 
E = cos A, the solution can be represented as 4>(x) n — 
P n {x) exp(— Xx)C(x), where P n (x) are some polynomials 
of degree n 
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In the extreme case, where m out = Qout, the spectrum 
degenerates (does not depend on n) and E = 1. The ab- 
sence of such an extreme state in the quantum spectrum 
is in complete agreement with the analogous result for the 
mass spectrum of a Reissner-Nordstrom black hole 
This result can be interpreted in terms of the postulated 
third law of thermodynamics for black holes, according to 
which the extreme state of a black hole is unattainable. 
At the quantum level, this means that the transitions 
(decays) of a black hole to the extreme state are impos- 



sible 



19l421j. Note that the formal solution of the wave 



equation (TT5)) in the extreme case is 



(x) = d(x) +C 2 (x)x, 



where C\ and C2 are periodic functions with a period of 
1. These periodic functions can be expanded into Fourier 



series to give 

00 



<K*> 



k=0 



Ck exp(—2-Kkx) + ix <^fc cxp(— 27rfcx), (20) 
fc=o 

where we made the inverse transformation x — > ix and 
took into account the boundary condition at infinity. All 
of the unknown coefficients c& and dk can be found using 
the boundary conditions (JTTJ) . 

The case where M 2 < Q 2 nt may be considered. In 
this case, a- lri can take on two values: <j- m = ±1. For 
radii p < (Qout — M 2 )/2m out < 2m out , the value of er in 
is negative and a[ n = — 1. It is easy to write out the 
wave equation and to solve it. The solution for the wave 
function will have the form <fi(x) = (—l) x (/)(x), where 
4>{x) is solution (fT4"|). The corresponding mass spectrum 
will not change. 

Let us now consider the more general case, where only 
Qin = 0. To find the general solution of Eq. pUj) . we will 
now make a different change of variable, namely, x — > ix. 
This transformation means the rotation through 7r/2 in 
the complex plane. Equation (110p will then be rewritten 
as 

2m- m Mi 



(x-1) + 



+ i- 
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M 2 



Qout 



<)>(x + l) 



(x) = 2Ecf>(x). (21) 



The latter equation can also be expressed in terms of 
Meixner polynomials [15I ] and its solution is 

P x Y{x + 1) 



</>n(x) 



C{x)- 

xA" 



/3 2x+2 "r( 7 + x) 
T(x + 1 - n) 



(22) 



where 



= E-^/E 2 - 1, 7 = i2m in Af, C(x) =C(x+l), (23) 

and C(x) is also a periodic function with a period of 1. 
This function can be similarly expanded into a Fourier 
series. Let us write out the first two polynomials: 

4>(x) n =i = lP 2 ( 1 + x)-x}0 x C(x), (24) 

4>(x)n=2 = X 4 (J + X + 1)(J + X) 

-2(3 2 x(j + x)+x(x-l)]C(x). (25) 
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If, as above, we set x = Xi and 2^+1 = X% + 1 and sum 
the wave functions with the weight 

r( 7 + a;) 



then the above polynomials will be orthogonal (lsj : 



(26) 



^ cj) n {xi)(j) m {xi)p(xi) = 5 nm d 2 n , (27) 



where 



2 _ n!r(n + 7) 



(28) 



" /3 2 »(i-/? 2 )^r( 7 ) 

Making the inverse change of variable x — > —ix we ulti- 
mately obtain the polynomials 

00 

<fi(x) n= i = [/3 2 (7 - ix) + ix]P~ lx V" c fc exp(-27rfca;), 

fe= — oc 

0(s)n=a = /3""[/3 4 (7-^ + l)(7-^)+2/3 2 zx(7-za;) 

OO 

+ 1)] ^2 c k exp(-2nkx). (29) 



fc=-c 



It is easy to show that at 7 = 0, i.e., at mi n = 0, solution 
([29} of the wave equation will be transformed into solu- 
tion ([20)1. The discrete mass spectrum is now specified 
by the equation 



z(M 2 - Q 2 out + 2m in M/3) = 2n^ - 1. (30) 

The mass spectrum was found to be imaginary, because 
Hamiltonian ((SJ) will no longer be Hermitian at m m 7^ 0. 
There exists one degenerate case where the spectrum 
is real: at E\ = 1 and M 2 - Q 2 out + 2m in M(3 = 0. 
From these two equations we obtain the conditions for 
the parameters of the problem under which the mass 
spectrum is degenerate and real: m out — mj„ + M, 

This limiting case repre- 



M = -m in + yfi 



Q 
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sents the transition from the discrete mass spectrum to 
the continuous one. To make the Hamiltonian Hermitian 
at mi n 7^ 0, it is necessary to make the change of opera- 
tors A(x)B(p) -> \[A(x)B{p) + B*(p*)A*(x)], which, in 
our case, corresponds to the change 



1 /. d \ 1 



1 ( d 
exp I % — 

x \ ox 



d \ 1 



exp l" Z fe x 



(31) 



After this change of operators, the wave function for the 
Hermitian Hamiltonian is 



(x + i) 



<P(X ' 

M 2 



ni;„M 



(x + i) <p{x — i) 
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' out 



(x) = 2E4>{x). 



(32) 



This equation is much more complex. Its solution and the 
energy spectrum must depend on two quantum numbers 
that define the mass spectrum of the inner, m m , and 
outer, TO out , black holes. We managed to find only an 
approximate solution of this equation in the case where 
the shell mass M was a small parameter of the problem, 
but the masses m ou t and mj n were not small and of the 
same order of magnitude. For simplicity, we also assume 
that Qout = 0. At large radii p, we see from Eq. ((D) 
that the difference of the black hole masses outside and 
inside the shell in the approximation under consideration 
is also a small quantity of the order of M. Consequently, 
the parameter of the problem E is not small. Making, as 
above, the change of variable x — ¥ ix in the wave equation 
(|32[) . we will obtain an equation with a linear accuracy 
in small parameter M 



0(X + 1) 

+im\r,M 



-4>{x-i) 

<f>(x + l) <t>{x-l) 



(33) 



x - 1 



2Ecj>{x) 



We will seek a solution of Eq. ([M]) in the form 4>{x) — 
4>o{x) + My(x). In the zeroth order in small parameter 
M, the wave equation will be reduced to a simple differ- 
ence equation 



4> {x + 1) + (j)v{x - 1) = 2E4> [x), 



whose solution is 



Mx) = C 1 (x)f3 x +C 2 (x)~p x , 



(34) 



(35) 



where the quantities (3 and $ are defined by Eqs. (fT5)) and 
([23|). respectively, while the functions C\{x) and C^ic) 
are periodic with a period of 1. 

In the first order of smallness, the function y(x) satis- 
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fies the equation 



y(x + 1) - 2Ey(x) + y(x - 1) = 



(36) 



Cx(x) 



X X—l 



' QX+1 p x ~ 1 



-C 2 (x) 



x x— 1 



whose solution is 



y{x) = 



^=A-C x {x) [^[2F(x, 1, 1+x, P 2 ) -1] 



+(3 X ~ 1 [f3 2 V(x + l) + <f{x)] 



+C 2 (x) \^L[2F(x, 1,1+x, /3 2 ) -1] 
L x 



+P x - i [p 2 ^(x + l) + y(x)] 



(37) 



Here F(x, 1, 1 + x, (3 2 ) and F(x, 1, 1 + x, (3 2 ) are the hy- 
pergeometric functions, while fy(x) = T (x)/T(x) is the 
logarithmic derivative of the gamma function. 

Consider the physical meaning of this solution. First, 
solution (|35|) represents two waves traveling in opposite 
directions. This can be easily made sure if the solution 
is represented as 



\I> ~ exp(-iEt ± x In (3), 



(38) 



which follows from the time-dependent Schrodinger equa- 
tion 



= Ety = ih- 



dt 



(39) 



The general solution (|35|) resembles the solution of the 
Schrodinger equation in a periodic potential (the Bloch 
function). However, this is only a qualitative similarity, 
because there are originally no periodic initial conditions 
in the problem under consideration. Since the solution 
obtained describes a continuous mass spectrum, the ap- 
proximate solution of the original equation (|34[) also has 



a continuous mass spectrum. We failed to find a solution 
with a discrete mass spectrum. As we see from the gen- 
eral solution (|3"7|) . it is also periodic with a period of 1. It 
can also be seen that the wave structure of the equation 
remains unchanged. In contrast, the amplitude is now 
a function of the argument and can change according to 
the presented law due to the addition of a nonlinearity 
to the difference equation. 



We considered here the simplest model of semiclassical 
quantization of a thin shell in the Reissner-Nordstrom 
metric. We used the hypothesis, based on the calcula- 
tions of the mass spectrum for a simpler model, that 
the mass spectrum of the shell must depend on two 
quantum numbers in the case where the metric is the 
Schwarzschild one inside and outside the shell. This can 
be explained in a simpler language as follows. Consider 
the Schwarzschild black hole. There exist two space-time 
regions R± on the Carter-Penrose diagram for this met- 
ric. In accordance with the results of p], 0], we will find 
for each of the regions R± that the black hole mass spec- 
trum must depend on two quantum numbers. 
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